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Abstract 

We show that conditional output measurement on a beam splitter may be 
used to produce photon-added states for a large class of signal-mode quan- 
tum states, such as thermal states, coherent states, squeezed states, displaced 
photon-number states, and coherent phase states. Combining a mode pre- 
pared in such a state and a mode prepared in a photon-number state, the state 
of the mode in one of the output channels of the beam splitter "collapses" 
to a photon-added state, provided that no photons are detected in the other 
output channel. We present analytical and numerical results, with special 
emphasis on photon-added coherent and squeezed vacuum states. In particu- 
lar, we show that adding photons to a squeezed vacuum yields superpositions 
of quantum states which show all the typical features of Schrodinger-cat-like 
states. 



1 Introduction 



It is well known that according to von Neumann's projection principle |l 
conditional measurement may be a fruitful method for quantum-state ma- 
nipulation and engineering. In particular, when a system, such as a correlated 
two-mode optical field or a correlated atom-field system, is prepared in an 
entangled state of two subsystems and a measurement is performed on one 
subsystem, then the quantum state of the other subsystem can be reduced 
to a new state. Systems that have typically been considered are the waves 
produced by parametric amplifiers 0, [3], |J and degenerate four-wave mixers 
M, the interfering fields in the output channels of a beam splitter M, [|, and 
systems of the Jaynes-Cummings type in cavity QED [|7| or trapped-ion stud- 
ies || . Further, state reduction via continuous measurement has also been 
considered || [10], [□], £|]. Conditional measurement offers new possibilities 
of generating extremely nonclassical states, such as photon-number states 
|, |, 0, 0, [TJ and Schrddinger-cat-like states |, 0, |T§. 

In this paper we show that zero-photon conditional output measurement 
on a beam splitter can be used advantageously to generate photon-added 
states for a large class of quantum states of the signal mode, such as ther- 
mal states, coherent states, squeezed states, displaced photon-number states, 
coherent phase states etc.. It can be expected that repeated application of 
the photon creation operator to the signal-mode quantum state can produce 
extremely nonclassical states. Photon adding can therefore be expected to 
improve the performance of noise reduction schemes In particular, the 
states that are obtained from coherent states by repeatedly applying to them 
the photon creation operator - the photon-added coherent states - can be 
regarded as non-Gaussian squeezed states and were first introduced and dis- 



cussed in |16| and it was proposed that they can be produced in nonlinear 
processes in cavities. We show that photon-added coherent states can also 
be generated via conditional output measurement on a beam splitter, mixing 
a signal mode prepared in a coherent state with a second mode prepared in 
a photon-number state. 

We further show that when a signal mode prepared in a squeezed vac- 
uum is mixed with photon-number states and zero-photon conditional output 
measurements are performed, then photon-added squeezed vacuum states 
can be produced. We analyze the states in terms of the photon-number and 
quadrature-component distributions and the Wigner and Husimi functions. 
We show that photon-added squeezed vacuum states exhibit all the typical 
properties of Schrodinger-cat-like states, so that photon adding can be re- 
garded as a method for producing Schrodinger cats. In particular, they are 
shown to be superpositions of two non-Gaussian squeezed coherent states 
that tend to Gaussian squeezed coherent states for sufficiently large num- 
ber of added photons. It is worth noting that although the scheme bears 



some resemblance to that in |14 |, the two schemes are quite different from 
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each other. That concerns both the second input quantum state (which is 
the vacuum in []14| ) and the conditional output measurement (detection of a 
nonzero number of photons in |nj). 

The paper is organized as follows. In Sec. ^| the basic scheme is explained 
and the conditional output states are derived. The possibility of the gener- 
ation of photon-added states is studied in Sec. |3|, with special emphasis on 
photon-added coherent and squeezed vacuum states. The problem of mixed 
photon-added states is addressed in Sec. f|. A summary and concluding re- 
marks are given in Sec. g. 



2 Basis equations 

It is well known that the input-output relations at a lossless beam splitter 
can be characterized by the SU(2) Lie algebra |T7|, T?|. In the Heisenberg 
picture, the photon destruction operators of the outgoing modes, bk (k = 1,2), 
can be obtained from those of the incoming modes, a^, as 

2 

b k = ^2 Tk,k' hk', (1) 

k'=l 

where 

(T \-Jn( cosfle^ sinfle*^ \ ( ) 

is a SU(2) matrix whose elements are given by the complex transmittance T 
and reflectance R of the beam splitter, 

T = cos 9 e iipT , R = sin 6 e iipR . (3) 

Using the Schrodinger picture, the photonic operators are left unchanged, 
but the density operator is transformed. In this case the output-state density 
operator g out can be related to the input-state density operator g- m as 

Qont = V^QtaV, (4) 

where V can be given by |17|, [18 



y — e -i(<PT-<PR)L 3 e -2i8L 2 e -i(<p T +ip R )L 3 ^ 

with 

L 2 = ^(«i02 - a\ai), L 3 = \(a\ai - a} 2 a 2 ). (6) 

Note that tp® is a global phase factor that may be omitted without loss of 
generality, y9 = 0. Applying elementary parameter-differentiation techniques 
|19|j , we can derive the operator identity 



g— 2iOL2 gtan(S)a2(ii g2 In cos(#)L3 ^— tan(6>)aja2 
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which [together with Eq. ([|)] enables us to rewrite V\ Eq. (||), as 



where h k = a[dk- 

Now, let us assume (Fig. |I|) that the modes that are fed into the first and 
second input channels of the beam splitter are prepared in a state described 
by a density operator £ inl and a Fock state g in2 = \n )(n \, respectively (for 
a review on the problem of the generation of Fock states, see |2(J). The 
input-state density operator can then be written as 

Qin(n ) = g inl <g> \n )(n \. (9) 

Using Eqs. (|J), (H), and (||), after some calculation the output-state density 
operator g ont = g O nt(n ) can be given by 

-i oo oo no no 

ft-S) = E E EE(«T 2+, ff 2+k 

\ 1 I n 2 =0 rn 2 =0 k=0 j=0 



\n 2 +m 2 



no\(n \fn -k+m 2 \fn -j+n 2 \ 

y/k\j\m 2 \n 2 \\\k J\j J\ m 2 )\ n 2 ) 

x T^df 2 (a\) k Q- ml d{{a\) n2 (T*)^ \n - k + m 2 ) (n - J + n 2 |. (10) 

It can easily be seen that when the second input channel is unused, i.e., n 
= 0, then Eq. ( |T0"|) reduces to the relation considered in ||, ||, [TJJ]. 

From Eq. (|10"D we see that the output modes are, in general, highly cor- 
related. When the photon number of the mode in the second output channel 
is measured and m 2 photons are detected, then the mode in the first output 
channel is prepared in a quantum state whose density operator f) uti(^cb m 2) 
reads as 

foutiCno, ma) = , , rr~ n. • (11) 

lr U \ m 2 1 Pout [ n ) \ m 2) ) 

The probability of such an event is given by 

oo no no 

P(n ,m 2 )=Tr 1 ((m 2 |eaut(T»o)|m 2 »= £ E E \R\^ +k ^\T\ 2 ^ + ^ 

ri\=\i—v j=/j, k=/j, 

(-iy +k n \nx\ fn -is\(n -is\(n 1 +j\(n 1 + k\ 

where 

v = n — m 2 , n = max(0, v). (13) 
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3 Generation of photon-added states 



Let us now assume that the (signal) mode in the first input channel is pre- 
pared in a mixed state 

0inl=£p*|$><$| ( 14 ) 

(J2$P<s> — 1, <p<s> < 1) and restrict attention to the events that no photons 
are recorded in the second output channel, i.e., 

m 2 = 0, v = n (15) 

in Eqs. (|TJj) - C[3~2f) . Note that such events can be detected using highly 
efficient avalanche photodiodes. Combining Eqs. (|T0| ) and (|11|) and using 
Eqs. ([14]) and ([15]), we find that the mode in the first output channel is 
prepared in a state 

Qouti(n ,m 2 = 0) = *no)(*nol , (16) 



l^n ) = ^= (a!)" T ftl |$), (17) 

/AC,., 



where 



'no 

M no being a normalization constant, 

N no = (^(T^a^ial^T^l®). (18) 

The probability of observing the conditional state £> uti(^o, m i = 0) can easily 
be found from Eq. (O) and reads 



P(n ) = P(n ,m 2 = 0)\R\ 2n ° £ |T| 2ni h + U °) (mlamlm). 



(19) 



The states |\1/ T10 ) , Eq. (|T7|), are obviously the conditional states observed 
in the case when the mode in the first input channel is prepared in a pure 
state |$). From Eq. ([17]) we see that for chosen n the conditional state 
l^no) is a photon-added state, with n photons being added to the state 
~T ni |$). In particular when the absolute value of the transmittance is close 
to unity, |T| ps 1, then the state |^) is - apart from a rotation in the phase 
space - close to the input state |$). To be more specific, let us consider the 
expansion of |$) in the Fock basis, 

oo 

l$>= E c «>i>> (2°) 

ni=0 

and assume that 

c ni for m > n max . (21) 
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We see that (apart from the rotation mentioned) ~ |$), provided that 



IT 



ni 



1 for rii < n n 



(22) 



Note that for any physical state |$) the expansion in Eq. (^) can always be 
approximated to any desired degree of accuracy by truncating it at n max if 
n max is suitably large. 

However, there are classes of states for which the corresponding photon- 
added states can be produced even when |T| is not close to unity. Let us 
consider a class of a parametrized states |$(a)) such that 



|$(«)) = ]T c ni (a)|m), 

ni=0 



with c ni (a) oc a 



in 



(23) 



Since in this case the relation |^) ~T ni |$(a;)) = |$(Ta:)) is valid, the state 
1^) obviously belongs to the class of states |$(o;)). Hence when the input 
state |$) belongs to the class of states |$(a)), then the conditional output 
state 

l*n ) ~ (a t ) no |^) ~ (a+H^Ta)) (24) 

is a photon-added states, with the photons being added to a state that 
also belongs to the class of states | $(«)). The extension of the above given 
considerations to mixed states is straightforward. Typical examples of classes 
of photon-added states that can be produced in this way are thermal states, 
coherent states, squeezed states, displaced Fock states, and coherent phase 
states. 

When the number of added photons, n , is sufficiently small compared to 
the photon numbers that mainly contribute to the state |^) and the photon- 
number distribution |(rii|\I/)| 2 is sufficiently slowly varying with m, then the 
photon-added state \^? no ) exhibits, in general, properties that are similar to 
those of the state Writing 

1/2 



ni+rio 



-"n\ 



■ T ni c 



ni 



m=0 C ni+n 

and using the approximations 

1/2 



c 



111 



c «i+no 

we approximately derive 



ri\\ 



n 



,n /2 



|ni+n ) (25) 



(26) 



(27) 



ni=0 

i.e., photon adding leaves the state |^) nearly unchanged. With increasing 
number of added photons qualitatively new properties are expected to be 
observed. To illustrate the method, let us consider photon-added coherent 
and squeezed vacuum states in more detail. 
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3.1 Coherent states 

Let us first assume that the input field is prepared in a coherent state, i.e., 
Qmi = |3>) ($|, where 



|*) = 



oo on 

-\/3\ 2 /2 P 

n=0 



(28) 



with /? = |/3|e^. According to Eq. the conditional output states then 
reads 



l*no> 



K 



where (3 1 = T(3 and 



ACo = n !L no (-|/3'| 2 ) 



(29) 



(30) 



L no (x) being the Laguerre polynomial. The states are photon-added coherent 
states and can be represented in the Fock basis as 



l*no) 



-|/3'| 2 /2 °° 



K 



E 



(n + n )\ 



1/2 



\n + n ). 



(31) 



From Eqs. (19) and fl28|), the probability of producing photon-added co- 
herent states is given by 



(32) 



Using the relations GT 



Ln (0) 



and 



X, L n 



n + n 



e 2 L nn (x - z) 



n=0 



n! 



(L" Q (x) is the associated Laguerre polynomial), we find that 
P(no) = |^" lW Ln.(-i^)- 



(33) 



(34) 



(35) 



In Fig. ^ the probability P{tiq) is plotted for two values of the beam-splitter 
transmittance |T|. If \T\ and/or hq are not too small, P{tiq) as a function 
of |/3 1 can attain a maximum, the position of which is determined by the 
(positive) solution of the equation 



|T| 2 L 



no \ 



l\2\ 



R\ 2 K,-A-W\ 2 )- 



(36) 



We see that even if or no are increased the probability -P(no) for detecting 
zero photons may increase due to destructive interference. In particular 
Fig. 0(b) reveals that the maximum is shifted towards larger values of \j3\ 
when n is increased. For example, assuming a highly transmitting beam 
splitter such that |T| 2 = 0.99, probabilities P(n ) of about of 10% may be 
possible, provided that and n are sufficiently large (\/3\ ~40, n ^20). 
Photon-added coherent states (|29D were introduced and studied in detail 
with special emphasis on their nonclassical properties, such a squeez- 
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m 

ing and sub-Poissonian photon statistics (for their quantum-phase statistics, 
see [0). In particular, in |]16[ analytical results for the Wigner and Husimi 
functions are given. We therefore may restrict attention to the quadrature- 
component distribution (i.e., the phase-parametrized field-strength distribu- 
tion) 

p no (x,(p) = \(x,ip\y no }\ 2 } (37) 

which can be measured in balanced homodyne detection. For this purpose 
we expand the eigenvectors \x, ip) of the quadrature components 



x(<p) 



-1/2 



-lip; 



in the photon-number basis |24 



exp 



1-2 
2 



£ 

n=0 



E n (x)\n) 



(38) 



(39) 



(H n is the Hermite polynomial). Using the identity |pH| 



k=0 



exp[2xz — z ) H n (x — z), 



(40) 



we find that 



y-n 



exp 



- v2\/3'\ cos((p + if pi 
H no (x-2- 1/2 |/3'|e i(v+ ^' ) 



(41) 



where tppi = <fp + <fT- Plots of p no (%, V 9 ) are given in Fig. ||for n = 1 [P(n ) w 
30%] and no = 4 [-P(no) ~0.84%]. The ir periodic narrowing (broadening) of 
the quadrature-component distribution reveals that photon-added coherent 
states can be regarded as some kind of squeezed states. Nevertheless, they 
are in general quite different from the two-photon coherent states [2S| widely 
used in squeezed- light description. It is worth noting that in contrast to the 
familiar two-photon coherent states the photon-added coherent states are 
non-Gaussian states. In particular, the variance of x((p) is given by 



(Ax(f)r = (x(ff) - (x(f)Y 
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2L„ (-2|/3f)2 



2|/5'r[L^ (-|/? / | 2 )L„, (-|/?'| 2 ) 

" Lnol-I/ 3 '! 2 ) 2 ] cos ( 2 (^ + <P?)) ~ 2|^i 2 ^ (-i^i 2 ) 2 

- L no (-|/3'| 2 ) 2 + 2(n + l)L no (-|/3'| 2 )L no+1 (-|/3'| 2 )}. (42) 

With increasing number of added photons, no, the squeezing effect is (for 
chosen \/3\' > 0) enhanced. Needless to say that when no = (i.e., when no 
photons are added), then the conditional state |^ ) simply reduces to the 
coherent state \/3'). Finally, it should be noted that when the vacuum is 
mixed with a Fock state and the zero-photon measurement in one of the 
output channels of the beam splitter is replaced with a measurement of the 
Q function (in perfect six- or eight-port balanced homodynings), then the 
conditional states of the other output channel are similar to the photon- 
added coherent states |J. 

3.2 Squeezed vacuum states 

Let us now consider a squeezed-vacuum input state 

|$) = |0) s = S(O|0), (43) 



where 



S(£)|0) = exp{-I - £*a 2 ] } | > = (1 - \k\ 2 ) 1/4 £ ^rj 1 « n |2n),(44) 



11/2 



n=0 



2™n! 



with £ = |£|e l</3 « and K = e tlpt tanh |£|. The conditional output states ( |i~T|) are 
seen to be photon-added squeezed vacuum states 

!*»«> 7=^( at )"°^(OIO>, (45) 

where f = \^'\ e i( -^ +2 ^ T \ and tanh = |T| 2 tanh In the photon-number 
basis they read as 

l^o) = 1 7=— £ W^') In), (46) 

where 

^ 1 r / n /1 ,\( n - n 0)/ 2 



W«0 = -TT7 V , 71 \ [1 + (-1)™] (47) 
T |(n - n ) + 1 

and k! = T 2 k. Note that when n = then |^n ) simply reduces to a squeezed 
vacuum state (fP|), with k' in place of k. 
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For the sake of tranceparency and without loss of generality we will as- 
sume that k' is real, i.e.,ip^ + 2ipT = kn, with k being an integer. Note that the 
effect of other phases is simply a rotation in phase space. Using the doubling 
formula for the Gamma-function, 

r(2,!) - 4 " r ( „), (48) 



r(n + l/2) 2v^F 
and the Gauss-series of the hypergeometric function 

V{a)Y{b)^ T(c + n) nV m 
the normalization constant M' no can be derived to be 



K = VT^no! F[i(n + 1), |(n + 2), 1; k' 2 \ . (50) 

From Eqs. (|IJ) and (fHj) [together with Eqs. (§£[) and (|9|)] the probability 
of producing photon-added squeezed vacuum states is derived to be 

P{n ) = |itf n °/L-H 2 F(n + 1, §, 1; k' 2 ) • (51) 

Examples of P(no) are plotted in Fig. |j. In particular we see that for not 
too small transmittance of the beam splitter (and chosen no) the probability 
P(no) can increase with the value of \k\, which is similar to the behavior 
shown in Fig. § for a coherent input state. Note that the mean photon 
number of the squeezed vacuum (4"1|) is given by |k| 2 /(1 — \k\ 2 )- 
Using Eq. (|£6| ), the photon- number distribution 

p no (n) = \(n\y no )\ 2 (52) 

of the photon-added squeezed vacuum states can be given by 

1 2 

— \bn,no(^')\ if n > n (53) 



and p no ( n ) — elsewhere, and 

AC = W F[i(n Q + 1), |(n + 2), 1; K ' 2 ] . (54) 

From Eqs. (^3]) and ([|7|) we easily see that when the number of the added 
photons, n , is even (odd), then the photon-number distribution is nonzero 
only for even (odd) photon numbers. This oscillating behavior of the photon- 
number distribution obviously reflects the fact that only even photon num- 
bers contribute to the squeezed vacuum state to which photons are added. 
In particular the mean photon number 

(h) =J2np no (n) (55) 
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may be rewritten as 



<n> = ^log 



~' no K 



(56) 



Using standard formulas for the derivative of the hypergeometric function 
[BEIl, we find that 



(h) = n + i« ,2 (n + l)(n + 2) 



i(n + 3) i(n + 4),2;/€- 



/2 



|(n Q + l) i(no + 2),l;K 



/2 



(57) 



As expected, for k' — > (n) approaches n and it increases with «' and n . 
Examples of the photon-number distribution are shown in Fig. |5] for n = 1 
[P(n )~23%] and n = 4 [P(n ) « 0.45%]. 

Let us now turn to the quadrature-component distribution. Combining 
Eqs. @7p, ©, and (|§), we derive 



where 



2 -n 



A^' v / 7rA ri o+ 1 



exp 



A 



a; 



H 



1 - n'e i2 v)/Ax 



A = 1 + k ,; 



2k' cos(2<^). 



(58) 
(59) 



Note that in the derivation of Eq. ( |58D the summation formula |2T 

(f + 4^)-"/ 2 - 1 / 2 



oo fc 

£ — H 2 fc +n (x) 
fc=0 K - 



x exp 



Azx 2 



H 



l + 4z/ n Vv / TT4i / 



z| < ~ 



(60) 



has been used. The quadrature-component distributions plotted in Fig. |5] 
correspond to the same parameters as in Fig. |5|. From inspection of Fig. ^ 
interference fringes for (p close to or it are seen, whereas for (p near 7r/2 two 
separated peaks are observed. The behavior is typical of a Schrodinger-cat- 
like superposition of two macroscopically distinguishable states. 
Next let us calculate the Wigner function 



j +oo 

W no (x,p) = - f dye 2ipy (x-y\^ no )(^n \x+y). 

71 J 



(61) 



Using Eqs. ( |39f) and (|46[) [together with Eq. ([601)1, after some calculation we 
obtain 



W no {x,p) 



2e 



tt^W"[2(k/ + 1)]"o 



+i 



x 



_x-y_ 



H 



»0 I 



(62) 
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where 



A 



1 + k' 

Performing the y integration |^TJ yields 

\K>\ n 0y/2 



(63) 



W no (x,p) 



7rAf^[2(l~^)]^+U2 exp { Xx2 A J 

no 
fc=0 



n Q 



(64) 



The Wigner functions in Fig. [7| are plotted for the same parameters as in 
Figs. H and [| We again recognize the typical features of Schrodinger-cat-like 
states. 

Finally, let us consider the Husimi function 

Q„ (:r,p) = ^|H* no >| 2 , (65) 

where \a) is a coherent state and a = 2~ 1 / 2 (x+ip). Recalling the expansion of 
the coherent states in the Fock basis, Eq. (p8|), and using Eq. (|46|) , we easily 
find that 



Qn (x,p) 



[x 



y 



2\n 



7r2 n o+W" 

J v no 



expj— | [(1 - k')x 2 + (1 + k')p< 



(66) 



Note that the Husimi function is a phase-space function that can be measured 
in multiport balanced homodyning, such as six-port [26| or eight-port p7 



detections. Since the Husimi function can be regarded as a smoothed Wigner 
function, the oscillating behavior and the negative values that are typical of 
the Wigner function (see Fig. |7]) cannot be observed. 

Schrodinger cat-like states are commonly defined as superpositions of two 
macroscopically distinguishable states. From Eqs. 
I^no) can be given by the superposition 



where 



with 



^o) = ^(l^: ) ) + l^; ) ) 
1 



and (|47|) it is seen that 
(67) 



I nn / 



AC (±) 



n=no 



n— no 



r[(n - n )/2 + 1] 

Using standard relations [ 2l| , the normalization factor AC (±) in Eq. 
be expressed in terms of derivatives of the hypergeometric function, 



• /v n 



& 



no 



Q K m 



-K 



{F(l,l,l; K ' 2 )+^'F(l,l 



3. 12 

'2' 



(68) 



(69) 



can 



(70) 
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and the normalization constant in Eq. flB"7|) is given by A = |(A/^'| ± - , /A/^ / ) 1//2 . 
To demonstrate the nonclassical properties of the component states 

in Fig. H we have plotted the Wigner function (x,p) of l^^) for various 
values of n . We see that with increasing n the Wigner function becomes 
less structurized and the negative values are more and more suppressed [e.g., 
for k' = 0.6 and n = 15 the Wigner function attains negative values of the 
order of magnitude of — 10~ 4 , Fig. |8|(c)]. It is worth noting that the compo- 
nent states l 1 !^) can be regarded as non-Gaussian squeezed coherent states 
that tend to the familiar Gaussian squeezed coherent states as no becomes 
sufficiently large. To illustrate the difference between the states and 
the Gaussian squeezed coherent states, let us consider the Husimi function 
{±)( nr n \ = | ^ a |v[/{±)^ |2^ w j^h a — 2~ 1 / 2 (x + ip). Using the expansions (^) 



and (F~ 



we derive 



\a\ 2n °e 



exp 



1 / ( 2 i *2\ 

±k la + a J 



Erfcf ^x/hda 



(71) 



Here, the relation [28 



E 



n= r(|n + l) 



exp(z 2 ) Erfc(— z) 



(72) 



has been used, where Erfc(z) is the complementary complex error function 
defined by 



2 t°° 

Erfc(z) = 1 - Erf (z) = — / dze' 1 

\/7l Jz 



(73) 



From Eq. (]71|) the asymptotic form of Q^> (x, p) for large no is derived to be 
(see |Appendix AT) 



n ! 



exp(— \a =F v^o| 2 ) ex P \ K ' ( ft2 + a * 



(74) 



From inspection of Eq. ([74]) we see that for large no the Husimi function 
Q^'(x,p) becomes a single-peaked Gaussian centered at iy^no, i.e., when 
n — >oo then the states tend to the familiar Gaussian squeezed coherent 

states. 



4 Mixed photon-added states 

So far we have assumed that a mode prepared in a photon-number state |no) 
is fed into one of the input ports of a beam splitter, so that exactly no photons 
can be added to the state of the (signal) mode fed into the other input port. 
In practice however, it may be more realistic to consider a statistical mixture 
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of photon-number states rather than a pure state, because of smoothings in 



the generation of photon- number states |2(J • It is worth noting that - apart 
from some smearing - the above given results remain valid as long as the 
statistical mixture of photon-number states is sub-Poissonian. Let us return 
to Eq. (|) and assume that 

ftnOo) = ftnl ® fln2, (75) 

where 

ftn2 = J2p^o K)(™o| ■ (76) 

no 

To be more specific, let us consider (as an example of a sub-Poissonian dis- 
tribution) a binomial probability distribution, 



Pn 



^p no (l-p) N - no if n <N (77) 



and p no = elsewhere (0<p<l). Note that for p^ 0, N — ► oo, and pN finite, 
the binomial distribution (|77|) reduces to a Poisson distribution, with n=pN 
being the mean photon number. Using Eqs. fl75| ) and (|76"|), in place of Eq. (|TJ|) 
we easily find that detection of no photons in one of the output channels of 
the beam splitter now yields the conditional (mixed photon-added) state 

£?outl(ra 2 = 0) = ^2pn Qoutl(n ,m 2 = 0) = ^ Pn P$\ *no)(*n | (78) 
"0 no,"!" 

in the other output channel. Accordingly, the probability of detecting the 
state is the average of P(n ) given in Eq. (|19|), 

P = J2Pn P(n ). (79) 

no 

Examples of the quadrature-component distributions p(x, ip) =J2n Pno Pn (x. 
of mixed photon-added coherent and squeezed vacuum states are plotted in 
Figs. H(a) and (b), respectively. In the figures it is assumed that p = 0.S and 
N = 5 [i.e., the mean photon number and the photon- number variance are n = 
4 and (An) 2 = 0.8, respectively]. Comparing Figs. ||(a) and (b) with Figs. |3|(b) 
and|6|(b), respectively, we see that the typical nonclassical features (such as 
squeezing and quantum interference) are preserved, even when the photon 
number state \n ) is replaced with a sub-Poissonian mixed state (|76|) (i.e., 
a smeared photon-number state). As expected, the probabilities ([79]) of ob- 
serving the mixed photon-added states become smaller than those obtained 
for pure Fock-state inputs [P ~ 0.07% for the mixed photon-added coher- 
ent state in Fig. ||(a) and P ~0.04% for the mixed photon-added squeezed 
vacuum state in Fig. |9](b)]. 
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5 Conclusion 



We have studied the problem of generating photon-added states using con- 
ditional output measurement on a beam splitter. When a single-mode radia- 
tion field is mixed with a mode prepared in a photon-number number state, 
then the mode in one of the output channels photon-added is prepared in a 
photon-added state, provided that in the other output channel no photons 
are detected. We have studied the conditions under which the photon-added 
states can be regarded as photon-added input states or states that belong 
the same class of states as the photon-added input states do. Typical exam- 
ples of input states to which photons can be added in this way are thermal 
states, coherent states, squeezed states, displaced photon-number states, and 
coherent phase states. 

Photon-added states are highly nonclassical states in general. In particu- 
lar, photon-added coherent states can be regarded as non-Gaussian squeezed 
states, as can be seen, e.g., from the derived expression for the quadrature- 
component distribution. Another interesting class of states that can be 
produced in the way described are photon-added squeezed vacuum states, 
which exhibit all the typical properties of Schrodinger-cat-like states. Pho- 
ton adding to a squeezed vacuum can therefore be regarded as a method for 
producing Schrodinger-cats. We have analyzed the photon-added squeezed 
vacuum states in terms of the photon-number and the quadrature-component 
distributions and phase-space functions, such as the Wigner and Husimi func- 
tions, and have presented both analytical and numerical results. We have 
shown that photon-added squeezed vacuum states may be regarded as su- 
perpositions of two non-Gaussian squeezed coherent states that tend to the 
familiar (Gaussian) squeezed coherent states as the number of added photons 
goes to infinity. 

With regard to possible experimental implementations, we have also per- 
formed calculations allowing for an input mode prepared in a statistical mix- 
ture of Fock states in place of a pure Fock state. As expected, mixtures of 
Fock states give rise to conditional states that can be regarded, in a sense, 
as smeared photon-added states. It is worth noting that when the photon- 
number distribution of the mixtures is typically sub-Poissonian, then the 
characteristic nonclassical features of photon-added states can be observed. 
In the paper we have demonstrated the effect of smearing assuming a bino- 
mial photon-number distribution. In particular, combining mixed Fock states 
of that type with a squeezed vacuum, the conditional output states also ex- 
hibit, in general, properties that are typically observed for Schrodinger-cat- 
like states. That is to say, the characteristic properties of squeezed vacuum 
states are (apart of some smearing) preserved. Clearly, when the photon- 
number distribution of the mixed Fock states becomes Poissonian, then the 
nonclassical properties that are typically associated with photon-added states 
disappear. 
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Appendix A Derivation of Eq. (|7g|) 



To find the asymtotoic form ( ff4| ) of Q^(x,p) defined in Eq. ([H]), we first 
consider the asymptotic behavior of |a| 2no e _ ' a ' 2 for no-^oo, 



I |2no e -|a| 2 ^ 



n l 
2im 

2miQ 



exp ^ 
exp 



\a 2 — n \ 2 \ 
2n ) 

\a + \/no\ 2 I" — \/no\ 2 



2n r 



(Al) 



which may be approximated by the sum of two Gaussians centered at a 



2m 



»(-|a-\/™o| 2 ) 



vx.})[ -|n - v ;/ | 2 ) + exp(-|a + A/n^| 2 



(A 2) 



Next let us consider the infinite-series approximation of the complementary 



complex error function [28 



Erfc(x + iy) ~ Erfc(x) + - — [1 — cos(2rq/) + i sin(2x?/)] 

2nx 



+ 



2g — x2 00 g — n 2 /4 



E 



[/n(z, y) + ig n {x, y)] , (A3) 



7T ~[ n 2 + Ax 2 
where f n (x,y) and g n {x,y) are given by 

f n (x,y) = 2x[l — cosh(m/) cos(2a;)] + nsinh(m/) sin(2a;|/) 

and 

g n (x, y) = 2x cosh(m/) sin(2x) + n sinh(ra/) cos(2xy). (A 4) 

From Eq. ( A 3| ) [together with Eqs. ( A 4 ) and ( |A 4 )1 we easily see that 

Erfc(x + iy) m Erfc(x) if \x\ > 1. (A 5) 

Further the approximations 



Erfc(x) s; if x > 1, 
Erfc(-x) « if x < -1. 



(A 6) 
(A 7) 
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are valid PB|| . Recalling Eq. 

1 



1|) corresponds to =F^ 



:K a. 



I]), we see that the argument x+iy in Eq. (|A| 
Since for chosen n f (V 7^ 0) and sufficiently large 
n the Gaussians in Eq. ([A 2|) are nonzero only for large \a\, we can apply 
the approximations (|A 5| ) - (|A 7[ ), which [together with Eq. ( |A 2|) ] yields the 
asymptotic form (|74l). Note that due to the error function only one of the 
two Gaussians in Eq. ( [A 2j ) survives. 
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Qinl 




£outl 



Figure 1: Scheme of the experimental setup. When a signal mode prepared 
in a state f) inl is mixed with another input mode prepared in a Fock state 
|n ) and in one of the output channels of the beam splitter no photons are 
recorded by a detector (D), then the quantum state f) out i of the mode in the 
other output channel "collapses" to a photon-added state. 
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(a) 




Figure 2: The probability of producing photon-added coherent states is 
shown as a function of \(5\ for two values of the beam-splitter transmittance 
[(a) |T| 2 = 0.3; (b) |T| 2 = 0.8] and various values of n . 
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(a) 




Figure 3: The quadrature-component distribution of photon-added coherent 
states for f3' = 0.89 = 1, |T| 2 = 0.8) and two values of n [ (a) n = 1; (b) 
n = 4]. 
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(a) 

0.8 | , 




Figure 4: The probability of producing photon-added squeezed vacuum states 
is shown as a function of \k\ for two values of the beam-splitter transmittance 
[(a) |T| 2 = 0.3; (b) |T| 2 = 0.8] and various values of n . 
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Figure 5: The photon-number distribution of photon-added squeezed vacuum 
states for k' = 0.6 (\k\ = 0.67, |T| 2 = 0.8) and two values of n (n = 1, full 
bars; n = 4, dashed bars). 
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Figure 6: The quadrature-component distribution of photon-added squeezed 
vacuum states for «/ = 0.6 (\k\ = 0.67, |T| 2 = 0.8) and two values of n [(a) n 
= 1; (b)n = 4]. 
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(a) 



W no {x,p) 




gure 7: The Wigner function of photon-added squeezed vacuum states for 
= 0.6 (\k\ = 0.67, |T| 2 = 0.8) and two values of n [(a) n = 1; (b) n = 4]. 
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Figure 9: The quadrature-component distributions of mixed photon-added 
coherent (a) and squeezed vacuum (b) states for p= 0.8 and N = 5 in Eq. (|77|) , 
the values of the other parameters being the same as in Figs. § and |(| 
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